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Hopf solitons in the Skyrme-Faddeev model are string-like topological solitons classified by the 
integer-valued Hopf charge. In this paper we introduce an approximate description of Hopf solitons 
in terms of elastic rods. The general form of the elastic rod energy is derived from the field theory 
energy and is found to be an extension of the classical Kirchhoff rod energy. Using a minimal 
extension of the Kirchhoff energy, it is shown that a simple elastic rod model can reproduce many of 
the qualitative features of Hopf solitons in the Skyrme-Faddeev model. Features that are captured 
by the model include the buckling of the charge three solution, the formation of links at charges five 
and six, and the minimal energy trefoil knot at charge seven. 



PACS numbers: 11.27. +d 

The Skyrme-Faddeev model Q is a field theory in 
three-dimensional space. It possesses finite energy topo- 
logical soliton solutions that have a novel string-like 
structure and are known as Hopf solitons, as the topo- 
logical classification is via the integer- valued Hopf charge 
Q. This is a linking number, rather than the more famil- 
iar degree that appears in most theories with topological 
solitons. 

It has been suggested Q that the model may be of rel- 
evance for a low energy description of QCD, where the 
solitonic strings could describe glueballs. Other applica- 
tions have also been proposed in the context of condensed 
matter physics, regarding the study of two charged con- 
densates [3j. 

Substantial numerical work [H-IH has produced a com- 
prehensive catalogue of minimal energy solitons in the 
Skyrme-Faddeev model. For Q = 1 and Q = 2 the mini- 
mal energy solitons are axially symmetric, but for larger 
values of Q they are more exotic and include knots and 
links. These solutions have been obtained using sophis- 
ticated, and computationally intensive, numerical simu- 
lations of the highly nonlinear field theory. The aim in 
this paper is to provide a simple approximate descrip- 
tion of Hopf solitons in terms of the string core, thereby 
reducing the three-dimensional field theory to an effec- 
tive one-dimensional description. The resulting theory is 
that of a generalized elastic rod, extending that of the 
classical Kirchhoff rod. It is shown that a minimal ex- 
tension provides a simple model that is already capable 
of reproducing many of the qualitative features of Hopf 
solitons. 

The Skyrme-Faddeev field is a map <fi : K. n- S 2 , 
which is realized as a three-component unit vector = 
(cpx,(j)2,4>3). As this paper is concerned only with static 
solutions then the model can be defined by its energy 



E SF = I ' di<t>- d t <j> + i(c^ x d.jct>) 2 d 3 x. 



(1) 



The first term in the energy is that of the usual 0(3) 



sigma model and the second is a Skyrme term, required to 
provide a balance under scaling and hence allow solitons 
with a finite non-zero size. 

Finite energy boundary conditions require that the 
field tends to a constant value at spatial infinity, which is 
chosen to be </>(oo) = (0,0,1). This boundary condition 
compactifies space to S 3 , so that the field becomes a map 
<f> ■ S 3 h-> S 2 . Such maps are classified by ^(S 2 ) = Z, so 
there is an integer- valued topological charge Q, the Hopf 
charge. It has a geometrical interpretation as the linking 
number of two closed curves obtained as the preimages 
of any two distinct points on the target two-sphere. A 
natural definition of the position of the soliton is pro- 
vided by the preimage curve of the point 4> = (0, 0,-1), 
which is antipodal to the vacuum value. The position of a 
Hopf soliton is therefore a closed string, or possibly a col- 
lection of closed strings since the preimage of any point 
may contain disconnected components. This is the novel 
string-like aspect of Hopf solitons. An energy bound of 
the form E > aQ 3 / 4 , has been proved [9j, where a is 
a constant, and it is known that the fractional power is 
optimal [10(. The sublinear growth has a simple physi- 
cal explanation, in that Hopf charge can be accrued by 
the linking of distinct components of the position string 
(or by self-linking of a single component), so that the to- 
tal Hopf charge can be greater than the naive sum of its 
components. 

Generically, the cross-section through the string posi- 
tion locally resembles a planar soliton of the 0(3) sigma 
model. The planar soliton has an internal phase, which 
may vary along the length of the curve to provide a twist. 
The simplest axially symmetric Hopf soliton of charge Q 
may be pictured as a circular string, with a constant 
twist rate, in which the total change in the phase angle 
around the circle is 2ttQ. Such solutions exist for all Q, 
but their energy grows linearly with large Q, in contrast 
to the Q 3 / 4 growth of the minimal energy solitons. Only 
for Q — 1 and Q — 2 are the minimal energy solitons of 
this axial form. For Q > 2 there is a buckling instability 
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and in fact the minimal energy soliton with Q = 3 has a 
buckled conformation. 

The Q — 4 minimal energy soliton is axially symmetric 
but is anomalous in that the cross-section is described by 
a double planar soliton, with total phase twist 4tt. There 
is a local minimum of the energy that describes a buckled 
Q = 4 solution but it has an energy that is a few percent 
above that of the global minimum. There is a further 
solution at Q — 4 that consists of two Q = 1 solitons 
that are linked once. If one links once a ring of charge Q\ 
and a ring of charge Q2, the resulting configuration has 
total charge Q = Qi + Q2 + 2 due to the extra linking of 
the preimage curves of the components. The charge four 
link has an energy between that of the buckled solution 
and the global minimum. 

Minimal energy solitons with Q = 5 and Q = 6 are 
both links, with two components and a single linking. In 
the first case the components have charges one and two, 
and in the second case both components have charge two. 
The first minimal energy knot appears at Q = 7. It takes 
the form of a trefoil knot, which has crossing number 
three, and has four units of twist, to make the total Hopf 
charge equal to seven. 

For later comparison, a selection of minimal energy 
Hopf solitons are displayed in the top row of Figure [1] 
For clarity the string position is displayed by plotting a 
tube around the string, given by an isosurface where </>3 = 
—0.8. The red curve indicates the twist and is plotted in 
a similar fashion from the preimage of a point close to 
the vacuum value. 





FIG. 1: Top row: The string position (yellow tubes) for min- 
imal energy Hopf solitons in the Skyrme-Faddeev model for 
charges Q — 1,3,6,7. The red curves indicate the twist. 
Bottom row: The corresponding minimal energy elastic rods 
(blue tubes). To aid visualization, the displayed rod thickness 
is only half the true thickness. 



The general form of an elastic rod energy will now be 
derived from the Skyrme-Faddeev held theory energy ([1} . 
The approach taken is to define a field configuration by 
specifying a closed curve 7(s) € R , where s € [0,L] 
is an arclength parameter, together with a quasiperiodic 
twist function a(s) such that a(L) = a(0) + 2nm, where 



m G Z. The field configuration is defined by introduc- 
ing tubular coordinates centered on the curve and then 
mapping the tube to the target two-sphere by first twist- 
ing each transverse disk using a(s), and then applying 
a fixed degree one map from the disk to S 2 . Physically, 
the curve may be pictured as the location of the string, 
the fixed degree one map is the analogue of the planar 
soliton and the twist function specifies its internal phase. 

To make the above construction explicit, let k(s) and 
r(s) denote the curvature and torsion of the curve ~f(s) 
and consider the standard Frenet frame with normal 
N(s) and binormal B(s). Take the radius of the trans- 
verse disk to be equal to the radius of curvature 1/k(s). 
The relation between the Euclidean coordinates x 6 M 3 
and the tubular coordinates s, p, 9 is 

x = 7(a) + ~(N(s) cos 9 + B(s) sin6>), (2) 

where p and 9 are polar coordinates in the unit disk. 
The degree one map from the disk to S 2 has the form 



4> = (sin / cos 9, sin / sin 9, cos /), 



(3) 



where f(p) is a profile function satisfying the boundary 
conditions /(0) = tt and /(l) = 0. 

Assuming the tubular region has no self-intersections, 
then for points in R 3 inside the tubular region the above 
construction provides an explicit field c/>, and outside this 
region the field is set to its vacuum value <j) — (0,0,1). 
This defines a field <f> throughout R with a well-defined 
Hopf charge. In the case that the curve 7 is an unknot 
then the Hopf charge is equal to the number of revolu- 
tions of the twist function, Q = m. 

Substituting this ansatz into the Skyrme-Faddeev en- 
ergy ([!]) gives, after a lengthy calculation, the energy 



E = 



bi + CiK 2 



C2 



^(&3 + C 3 K 2 ) }ds, (4) 



where the constants are given in terms of integrals of the 
profile function and its derivative. The terms in the en- 
ergy associated with the coefficients b\, 62, 63 derive from 
the sigma model term and the terms associated with the 
coefficients c±, C2, C3 are obtained from the Skyrme-term. 

The energy ^ is a function of the position curve ~f(s) 
and the twist function a(s), and is a natural generaliza- 
tion of the classical Kirchhoff energy for elastic rods. The 
Kirchhoff energy is recovered if all coefficients except C\ 
and C2 are set to zero. In this case, the two terms cor- 
respond to a bending and a twisting energy and the rod 
length L needs to be fixed. 

The main aim of the present paper is to show that 
an elastic rod model can reproduce many of the qualita- 
tive features of Hopf solitons. This motivates a limited 



3 



phenomenological approach, in which a simplified form 
of the energy ((4]) is considered, by setting some of the 
coefficients to zero. The remaining coefficients are then 
fixed by fitting to selected properties of Hopf solitons. 

The starting point is the Kirchhoff energy, therefore Ci 
and C2 are taken to be non-zero. To balance the scaling 
properties of these two terms requires at least one contri- 
bution that is obtained from the sigma model term. We 
choose the simplest possibility, namely to take b\ to be 
non-zero, so that there is a contribution to the rod energy 
that simply involves its length L. The reduced elastic rod 
energy is therefore defined by setting 62 = &3 = C3 = 0. In 
particular, this means that terms involving the derivative 
of the curvature are ignored in the reduced theory. 

The reduced elastic rod energy contains only three pa- 
rameters. Two of these parameters simply define the 
energy and length units of the model, which can be fixed 
by matching to the energy and length of the Q = 1 Hopf 
soliton. This means that, in suitable units, two of the 
three coefficients may be set to unity and the energy of 
the reduced elastic rod model can be taken to be 

E = J o (l + k 2 + C(a' -t) 2 ^ ds. (5) 

A rod that is a minimizer of the energy ([5]) is in fact a 
Kirchhoff rod, as it is also a minimizer of the Kirchhoff 
energy. The additional feature induced by the first term 
in |5]) is that the rod length L is not fixed, but is itself 
determined by energy minimization. 

Axially symmetric Hopf solitons of charge Q are mod- 
elled by circular rods with a linear twist function a = 
2nQs/L. They have energies and lengths given by 

E° = An^l + CQ 2 , L§ - 2ny/l + CQ 2 . (6) 

Note that this energy formula captures the linear growth 
for large Q found for axially symmetric Hopf solitons. 

As shown by Michell [1JJ in 1889 , a circular rod will 
buckle if the total twist 2ttQ exceeds a critical value 
2nQ > 27rv / 3/C. For the first buckling instability of an 
axial rod to occur at Q = 3 requires that -j= < C < 
The analysis below reveals that the rod energies under- 
estimate the soliton energies, so twisting is not suffi- 
ciently penalized. To minimize this deficiency C should 
be taken to be close to its upper limit, therefore in the 
following we set C — 0.85. With this value of G the 
above formulae yield Eg = 17.09, Lg = 8.55, Eg = 
26.36, Lg = 13.18, Eg = 36.96, Lg = 18.48. A com- 
parison of the Q — 1 and Q = 2 rods reveals that 
Eg I Eg = Lg/Lg = 1.54. For Hopf solitons these ratios 
are E^/Ef* = 1.63 and Lf/Lf F = 1.45, so the rod 
model is reasonably accurate, with an error of around 
6% for both the energy and length of the Q = 2 solution. 
The Q = 1 axial rod is displayed as the first image in the 
bottom row of Figure [TJ 



To avoid self-intersection of the rod we impose an 
additional constraint by assigning a radius, R, to the 
cross-section of the rod and demand that this thick- 
ened rod does not self-intersect. For the rod thickness 
to have no influence on the axial solutions requires that 
R < Li/(2n) = 1.36. To maximize rod energies, we set 
R equal to this upper limit. 

To numerically compute minimal energy rods the curve 
7 is discretized into a polygonal curve with 100 vertices. 
A simulated annealing method is used to minimise a dis- 
crete version of the energy (O , with a condition of equal 
edge lengths imposed using a penalty function. A discrete 
version of curvature is calculated from the angle between 
two neighbouring edges [13] and the self-intersection of 
the rod is excluded by imposing an upper limit on the 
curvature n < 1/R, and a lower limit on the separation 
between specific subsets of pairs of vertices, using the al- 
gorithm presented in [l2T| . The twist contribution to the 
energy is evaluated by identifying the combination a' — t 
with the rate of change of the angle between the mate- 
rial frame of the rod and the twist free Bishop frame, as 
described in [l3j |. 

Our numerical computations determine the energies 
of the first two buckled rods to be E^ = 35.35, and 
E™ — 44.16. The buckled minimal energy Q = 3 rod is 
displayed as the second image in the bottom row of Fig- 
ure [TJ where it can be seen that it closely resembles the 
corresponding Hopf soliton. 

Links can be formed by putting together two rods, pro- 
vided their radii are first increased to allow a sufficient 
gap to accommodate the other component of the link. 
Let Eq denote the energy of the axial rod with a radius 
at least 2R, so that there is enough room to slip another 
rod through its centre. Note Lg > 2Lg , hence for Q > 3 
the axial rod is already large enough to insert another 
rod giving Eq = Eq. However, for Q = 1 and Q — 2 the 
natural radius must be increased to give 

E$ = ^-(AR 2 + 1 + CQ 2 ). (7) 
^ R 

With the values of C and R given earlier this formula 
results in the energies Eg = 21.36, Eg = 27.25. Form- 
ing a single link between a rod of charge Q\ and a rod of 
charge Q 2 produces a link of charge Q = Q1 + Q2 + 2. We 
denote the energy of such a linked rod configuration by 
Eq 1 '® 2 . If both Qi and Q2 are less than three then both 
components of the link will be axial rods and the energy 
is simply obtained by addition as E^ 1 '® 2 — Eq i + E^ 2 . 

This gives the link energies E{' 1 = 2Eg = 42.73, E^ 1 = 
Eg + Eg = 48.62, E 2 ^ 2 = 2Eg = 54.51. Thus for Q = 4 
the link is slightly lower in energy than the buckled rod, 
which agrees with the result for Hopf solitons, where so- 
lutions of both types also exist. Recall that the minimal 
energy Hopf soliton for Q = 4 is anomalous, in that its 
cross-section is a double planar soliton, therefore in its 
current form the elastic rod model is unable to describe 
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this soliton. It is possible to generalize the rod model to 
a cross-section with arbitrary degree, but this refinement 
is more appropriate for future investigations that include 
all the possible terms in the energy. A lower bound on 
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FIG. 2: The data points are (normalized) elastic rod energies: 
circles for circular rods, triangles for buckled rods, squares for 
links, and diamonds for knots. The upper solid curve is the 
(normalized) minimal energy of Hopf solitons in the Skyrme- 
Faddeev model. 

the energy of a link formed from rods of charge three and 
two is Ej' 2 >E%>+E$ = 62.60. The energies of trefoil 
knots with charge Q are denoted by Eq , and have Q — 3 
turns, with a contribution of three to the Hopf charge 
due to the self-linking of the trefoil knot. Our numerical 
computations provide the following trefoil knot energies 
Ef = 53.09, Ef = 57.20, E? = 61.32. 

Comparing the above link and knot energies shows that 
links are preferred for Q = 5 and Q = 6, but at Q = 7 the 
knot energy is below the lower bound of the link. These 
results show that the forms of minimal energy rods for 
Q = 5,6,7 are in agreement with those of minimal energy 
Hopf solitons. The final two plots displayed in the bottom 
row of Figure Q] are the minimal energy Q = 6 and Q = 
7 rods, whose forms compare reasonably well with the 
corresponding minimal energy Hopf solitons presented in 
the top row. 

The data points in Figure [2] are the various rod ener- 
gies, normalized by the energy Ef of the Q = 1 rod. For 
comparison, the upper solid curve is the minimal energy 
in the Skyrme-Faddeev model, normalized by the energy 
Ef F of the Q — 1 Hopf soliton. Although the rod energies 
grow too slowly in comparison to the soliton energies, we 
have seen that several qualitative features are correctly 
reproduced regarding the minimal energy configurations, 
that is, axial for Q — 1,2, buckled for Q — 3, links for 
Q = 5, 6 and the first minimal knot at Q = 7. 

It is perhaps not surprising that the rod energies under- 
estimate the Hopf soliton energies, because the non-local 
interaction between different parts of the rod is very sim- 
plistic and does not provide a direct contribution to the 



energy. It might be possible to produce a more accurate 
model by introducing a more sophisticated non-local in- 
teraction, but then some of the simplicity and elegance 
of the current model would be lost. Similarly, the ex- 
tra terms neglected to produce the reduced rod model 
(for example, terms involving the derivative of the cur- 
vature) could also be considered and will again lead to 
increases in rod energies. This may provide a more ac- 
curate quantitative rod model, and could be investigated 
in the future. 

It is certainly of interest to obtain an improved quan- 
titative description of Hopf solitons, to build on the suc- 
cessful qualitative results presented here. An accurate 
rod approximation would be a very useful tool for inves- 
tigating Hopf soliton issues that are currently difficult 
to study within the field theory, such as the existence 
of non-torus knots and the structure of solitons for large 
Hopf charges. 

Finally, note that if the rod energy simply consisted 
of the first term in ([5]) then our problem would coincide 
with the construction of ideal knots and links [14| . in 
which the energy function is the length of the rod for 
a fixed thickness. The energy ([5]) is therefore an inter- 
esting hybrid of the Kirchhoff elastic rod and ideal knot 
energies. 
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